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Modeling active galactic nuclei: ongoing problems for the faint-end 
of the luminosity function 
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ABSTRACT 

We consider simple semi-analytic models that relate the active galactic nuclei (AGN) 
evolution to the merging history of their host dark matter haloes and quantify their ability 
of matching the AGN luminosity function and its spatial clustering at low and intermediate 
redshifts. In particular, we focus on the recent determinations of the AGN luminosity function 
in the hard X-ray band at z ~ which constitutes the most stringent observational test for 
our models. Indeed, while we find an acceptable agreement between the model bolometric 
luminosity function and the data at 1 < z ^ 2 and for luminosities larger than 10 10 Lboi,©> no 
semi-analytic model is capable of reproducing the number density of faint X-ray sources in 
the local universe. Some improvement can be obtained by advocating energy feedback that 
we model through a time-dependent Eddington ratio. Even in this case, however, the number 
density of faint AGNs is significantly below observations. This failure indicates that major 
mergers cannot constitute the only trigger to accretion episodes in the local AGN population. 
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1 INTRODUCTION 

The luminosity function (LF) of AGNs, namely the deriva- 
tive of their co-moving number density with respect to lumi- 
nosity, is a very important tool to understand their cosmologi- 
cal evolution. In the recent years, the AGN LF has been mea- 
sured in a wide range of redshifts, lu minosities and in differ- 
ent wavelen gth bands: radio (see e.g. iNagar et"aL 20051), opti- 
cal (s e e e.g. iKoehler et al.ll 19971: iGrazian et al.ll200C : Croom et al.l 
120051 : iRichards et al.1 120051 120061: ISiana et alj l2006h . infra-red 
(see e.g. iBrown et al.l 120061; iMatute et alj|2006l) . soft X-ray (see 
e.g. jMivaji et al. , uOOll: lHasinger et al.l l2005h. hard X-ray (see 



Ueda et al.l2003: Sazonov & Revnivtsev 2004l: lLa Franca et al.l 



120051 : IShinozaki et alj|200fj : iBeckmann et alj|2006h . In this work 
we are mainly interested in the two most recent determinations of 
the AGN LF in the hard (^ 2keV) X-ray band that, despite be- 
ing very local, provi des strong constraints to AGN models. The 
first one, provided by Shinozaki et al. (2006) (hereafter S06), con- 
sists of a complete, flux-limited sample of 49 sources from the 
HEAO-1 All-Sky catalogue, complemented with spectral informa- 
tion from ASCA, XMM-Newton and Beppo-SAX observations. 
All objects in the catalogue are optically classified as emission-line 
Seyfert galaxies at high galactic latitude (b ^ 20°) with column 
density N H > 10 215 crrr 2 and L x = L[2 - lOkeV] > lO^ergs" 1 . 
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The second AG N LF has been de t ermin ed in a harder X-ray band 
[20 - 40keV] by IBeckmann et al.1 J2006I) (hereafter B06), using a 
sample of 38 objects, preferentially located at low galactic lati- 
tude, detected by the imager IBIS/ISGRI on-board INTEGRAL, 
with Lx = L[20 — 40keV] > 10 41 erg s~ 1 . The main reason for con- 
centrating on these two datasets is that they allow to span a large 
band [2 — 40keV] in which it is possible to detect absorbed and un- 
absorbed AGNs at all galactic latitudes. 

Further observational constraints to theoretical models are 
provided by the AGN spatial clustering, which is often quantified 
by means of the angular or spatial two-point correlation function. 
Uncertainties in current modeling of the AGN clustering, however, 
make this second constraint less effective than the LFs. In spite of 
that, in this work we also check the ability of our models in match- 
ing the AGN biasing function, defined as the ratio between the spa- 
tial two-point correlation function of AGNs and dark matter (DM): 
b 2 (M,z) = £,agn I^dm, where M is the mass of the hosting halo. 
Also in this case, we consider the very recent determination of the 
AGN biasing function in the B-band bv lPorciani & Norberd fcOOfJ) 
(hereafter PN06), which provide the most accurate estimate of the 
AGN clustering to date. 

Ove r the years, sever al analyti c ( see e.g. lEfstathiou & Reesl 
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logical constant, Q.\ = 0.7, and a CDM density power spectrum 
with primordial spectral index n = 1 . 



_ and semi- 
analy tic ( see e.g. ICattaneo et al .1 119991; iKauffmann & Haehneltl 
20001; ICavaliere & Vittorinil I2OO0I; ICattaneol I200T 
Cavaliere & Vittorini 2002; Enoki et al. 2003; Volonteri et al 
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2006) oposed to describe 

the co-evolution of supermassive black holes (BH) powering 
AGNs and their DM halo hosts within the hierarchical clustering 
framework. These models proved to be successful in matching 
the optical LF of AGNs around their peak of activity at z ~ 2. 
However, modifications to the simple analytic models in which 
the AGN activity is triggered by haloes' mergers have been 
subsequently introduced to bring predictions into agreement with 
the observed AGN LFs both at higher and lower redshifts. For 
example, selective B H accretion at early stag es occurring either at 
super-critical rate (e.g Volonteri & Rees 200a) or at Eddington rate 
for B Hs hosted in the high-density peaks (e.g. IVolonteri & Reesl 
2006), can explain the observed number density of AGNs at z — 6. 

At low redshifts (0.5 < z ^ 2), inefficient cooling in large 
haloes is required to improve the match to the bright end of the LF 
in both the optical and hard X-ray bands (Marulli et al. 2006) (here- 
after M06). In the same redshift range. IVolonteri et al.l d2006h have 
recently shown that inefficient accretion, suggested by the results 
of numerical simulations that track accretion onto BHs following 
halo mergers ( Hop kins et ai]|2005t) . increases the number density 
of faint AGNs matching their optical luminosity function. Whether 
these modifications can also reproduce the LF of AGNs in the very 
nearby universe and in other luminosity bands it is still matter of 
debate. 

Before addressing this question, it is worth stressing that com- 
parisons between model and observed LFs are hampered by the 
fact that, while the former refer to bolometric luminosities, the lat- 
ter are measured in some specific bands. Luminosity band correc- 
tions represent therefor e a key issue tha t has been addressed by 
several authors (see e.g. lElvis et al ] |l994l ; lMarconi et alj[2004h . In 
thi s paper, we adopt the bolometric corr ection recently proposed 
by Hopkins. Richards & Hernquist (2006) (hereafter H06) using a 
number of determinations of AGN LFs in the interval < z ^ 5 
and in different bands ranging from infrared to optical, soft and 
hard X-ray bands. It is worth noticing that the H06 bolometric cor- 
rection was not calibrated using the S06 and B06 results. Yet, as 
we have verified, the bolometric correction successfully applies to 
both LFs. This is less obvious for the B06 than for the S06 dataset 
which, in fact, contains a significant fraction of objects that have 
been included in the H06 analysis. 

The main aim of this work is to predict the LF of AGNs using 
standard semi-analytic hierarchical models at z ~ and compare it 
to the S06 and B06 results by using the H06 bolometric corrections. 
In doing so, we also show model predictions for the AGN LFs and 
clustering up to z = 2, hence updating the results of M06. 

The outline of the paper is as follows. In the first section we 
briefly present the semi-analytic models considered in this work 
and summarize the main assumptions used therein. In Section 3 
we compare model predictions with the observed luminosity and 
biasing function of AGNs. Finally, in the last section we discuss 
our results and draw our main conclusions. 

Throughout this paper we assume a flat ACDM cosmologi- 
cal model with Hubble constant h = /fo/lOOkms -1 Mpc -1 = 0.7, 
a dominant contribution to the density parameter from the cosmo- 



2 SEMI-ANALYTIC MODELS 

Most of the semi-analytic models used in this work have been al- 
ready considered by M06. They all assume that the evolution of 
AGNs can be followed within the standard hierarchical scenario 
of structure evolution under the hypothesis that, at every redshift, 
the AGN activity is solely determined by the cosmological merg- 
ing history of their DM halo hosts. In this framework, the merg- 
ing history of dark haloes is described by the extended P ress & 
Schechter formalism dBond et alJl99ll ; lLacev & Colell993l) . while 
phenomenological prescriptions are adopted to model the feeding 
of the central BH and the physical processes of the AGN activity. 

Like in M0 6, the models are based on the s emi-analytic code 
developed by IVolonteri. Haardt & MadaJ fc003l) . AGN activity is 
assumed to be triggered when major mergers occur, i.e. when the 
two merging haloes have a mass ratio larger than 0.1. Due to the 
lack of an exhaustive study of the ultimate consequences of a 
galaxy merger in its whole parameter space, we are forced to make 
some simplifying assumptions to follow the merging events. Fol- 
lowing Cox (2004), we can assume that all halo mergers, except 
the ones with mass ratio smaller than 0.1, can destabilize the gas 
at the centre of the more massive halo, and consequently induce 
star formation and BH mass accretion. Notice that this threshold 
is lower than the value of 0.3 generally used in the literature, but 
it is not low enough to reproduce the observed faint AGN number 
counts, as we will describe later. So, an higher value of the mass 
ratio would worsen our resu lts, while a lower one would be in dis- 
agreement with the results o flCoxl ( l2004l) . which show that a typical 
merger with a mas s ratio of 0.05 does not induce starbursts. More- 
over, according to lTaffoni et alj ( 120031) , when P < 0.1 the dynami- 
cal friction timescale is larger than the Hubble time, hence prevent- 
ing the merging of satellite galaxies and, reasonably, making the 
accretion efficiency onto the central BH very low. 

We stress here that considering galaxy rather than halo merg- 
ers could lead to different predictions as it is not always true 
that a major merger of two dark matter haloes will result in a 
major merger of their galaxies (and viceversa). A galaxy-merger- 
driven scenario has the virtue of matchin g the observed correla- 
tion between BHs and their galaxy hosts dMagorrian et alj 1 19981 ; 
iFerrarese & Merrittl 12000 ; Geb hardt et ail ! 2000) but can only be 
self-consistently implemented within the framework of a full semi- 
analytic or numerical model of galaxy formation and evolution. 
This is beyond the scope of this work in which, instead, we use 
a model meant to minimise the number of free parameters. 

It is also assumed that the seed BHs formed with masses 
of 15OM0 following the collapse of the very rare Pop III stars, 
in minihaloes forming at z = 20 from th e density peak s abov e 
a 3.5o threshold; however, as shown by Volo nteri et alj (2003), 
the final results are not very sensitive to this choice. After 
every major merger, the BH at the centre of the more mas- 
sive halo increases its mass, Mbh, after a dynamical free-fall 
time when a significant fraction of the gas falls to the cen- 
tre of the merged system |Springel.~D i Matteo & Hernquistj |2005l ; 
|Pi MatteoTS pringel & Hernquist 120051) and is accreted at an ap- 
propriate rate. To implement this mechanism we need to specify 
the prescriptions for the mass accretion, that can be modeled by the 
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following relation 

M BH (f + 8f) =M BH (f)expn 



8t 1 - £ 

—fEdi(t)—— 
'Edd e 



(1) 



in which r^dd = 0.45 Gyr, /Edd is the Eddington ratio of mass ac- 
cretion, and e is the radiative efficiency. 

We have considered four different scenarios. 

• In the first one, labeled El in all plots, the accreted mass is 
proportional to the mass of the available gas and hence to the to- 
tal mass of the more massive progenitor: AM accl - = (xMh a i . We set 
a = 2 x 10~ 5 , in agreement with the normalization of the Mbh — 
relation at z = 0, where O p is the velocity dispersion of the host 
galaxy dTremaine et aljl2002l). scaling w ith the halo circular veloc- 
ity, v c , as suggested by iFerraresd d2002h . As M(j a i °= v^, the slope 
of the Mbh ~ rela tion is flatter than the observed one (but see 
IWvithe & Loebll2005h . 

• The second scheme, labeled E2, assumes a scaling relation 
between the accreted mass and the circular velocity of the host 
halo, AM accr oc fe-vj, which is normalized a-posteriori to repro- 
duce the observed relation between Mbh an d v c at z = ( Ferrarese 
2002). As in M06, we assume a linear dependence of k on red- 
shift, as k(z) = 0.15(1 + z) +0.05, in order to account for the de- 
crease of the gas available to fuel BHs. Unlike model El, here 
and in the other three models the relation between Mbh an d Mj, a i 
evolves in redshift as in IWvithe & Loebl J2003l) : M BH <* Mj^, ■ 
(1 +zf /2 ■ (A c /fl,„(z)) 5 / 6 , in which A c (z) = 18rc 2 + 82d - 39d 2 , 
d = £l m (z) — 1 and £2 m (z) represents the mass density parame- 
ter. Finally, we account for inefficient cooling in large haloes by 
preventing accretion within haloes of masses M na i > 1O 13 ' 5 M . 
It is worth noticing that this prescprition has a physical motiva- 
tion connected to both galaxy and AGNs formation since it has the 
same effect of including the low luminosity radio mode AGN heat- 
ing, as done in many semi-analytic models of galaxy formation to 
produce a massive galaxy population similar to the one observed 
(see e.g.lKang et alfooft iBower et alj|2006l : ICattaneo et alll200d : 
ICroton et alj2006h . 

• Another model, labeled B, assumes an early stage of super- 
critical accretion during which the central BH (AGN) accretes 
mass at a rate that can be estimated by the Bondi-Hoyle formula 
(Bond i&Hoyielll944ri . This model applies to metal-free haloes, 
therefore we assume that by z = 12 the interstellar medium has 
been enriched, and we inhibit super-critical accretion rates. When 
the super-critical phase ends, accretion proceeds in subsequent 
episodes as in model E2. This po ssibility has been recently advo- 
cated by Volonte ri & Reesl ( feOOSl) to reconcile a hierarchical evolu- 
tion with the existence of QSO at z ~ 6, hosting SBHs with masses 
~ 10 9 M @ . 

In all these models, which do not account for possible feedback 
mechanisms, we set in Eq. (TJ the Eddington ratio of mass accre- 
tion, / Edd = M/M Edd = 1. 

• In this work we consider also a different model, labeled H, 
which accounts for the results of the recent hydrodynamic sim- 
ulations of galactic m ergers in which AGN feedback is included 
( Hopk ins et alj 120051) that show that the Eddington ratio is not 
constant but depends on AGN luminosity. As the main variable 
in the model is the BH mass rather than the AGN luminosity, 
I Volonteri et al.| j2006) model /E dd (/) as follows. First the time spent 
by a given AGN per logarithmic luminosity interval is approxi- 
mated as 



where t Q ~ 10V, a = -O.95 + O.321og(L peak /lO 12 L ), and Lp eak 
is the luminosity of the AGN at the peak of its activity and 
it can be approximated as the Eddington luminosity of the BH 
at its final mass, i .e. when it sets on the «j B h ~ Gg relation 
(Hop kins et"ai]2 005). Then, since the AGN luminosity can be writ- 
ten as L = £/£d£/(f)ME dd c 2 , where £ is the radiative efficiency, 
£ = L/(/Edd( ? )MEdd c ' 2 ) = O.lQ, the following differential equation 
is used to describe the evolution of /E dd (?) 



d/EddO) _ / E dd a (0 ( EMEddC 



dt 



|a|f fi V 10% 



(3) 



dlogL 



\a\t Q 



10 9 L, 



(2) 



The instantaneous Eddington ratio /EddW ls obtained by solving 
Eq. [3] Model H assumes that the final mass of the black hole is 
determined by the circular velocity of the host halo, as in model 
E2. 



3 MODELS VS. OBSERVATIONS 

3.1 The bolometric luminosity function between 0.5 < z < 2 

To compute the model LF we have considered binary merger 
trees with masses in the range (1.43 x lO n M , 1O 15 M ) 
dVolonteri et~al]|2003l) . M06 have used models El, E2 and B de- 
scribed in the previous section to follow the accretion history of the 
BHs and the evolution of AGN luminosity at the centre of haloes. 
Here we have repeated this procedure to implement the novel model 
H . 

We have calculated the bolometric LF by simply discretizing 
the luminosity range of our modeled AGN sample in finite bins, 
counting their number density in all our merger tr ees, each of them 
weigh ted by the halo number density at z — ( Shet h & Tormenl 
1999) and normalizing for the number of merger trees and time- 
steps considered. Uncertainties in the model LF have been com- 
puted assuming Poisson statistics. 

As a first step, we extend the M06 analysis by comparing the 
bolometric luminosity function of H06 with the analogous quantity 
predicted by the four models at redshifts z = 0.5, 1.0, 1.5 and 2. 
This comparison represents a more severe test to the models than 
the one performed by M06, since the number of AGNs used by 
H06 to model their bolometric LF is significantly larger than those 
considered by M06, and consequently the uncertainties are smaller. 
The results are shown in Fig. Q] which is divided in four sets of 
plots, each one referring to a different model. In each set, com- 
posed by four panels, the dark-shaded area represents the la un- 
certainty strip around the model LF. The dashed curve shows the 
bolometric LF of H06 along with la uncertainty strip, plotted as 
a light-shaded area. In all plots the vertical, dashed line shows the 
minimum bolometric luminosity accessible to observations, £<J!u?, 
which turns out to be remarkably constant in the interval of red- 
shifts considered. Model predictions extend up to a maximum lu- 
minosity ^™ de i ~ 1O 14 L resulting from having set an upper limit 
to the mass of our DM haloes (M max ~ 1O 15 M ). 

To quantify the consistency between models and data we have 



' The radiative efficiency has been self-consistently determ ined by track- 
ing th e evolution of BH spins throughout the calculations iVolon teri et alj 
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Table 1. Values of S model and the corresponding per cent probability p(s zamode '), indicated in parentheses, for each models and redshifts. 



Model 






E(P(S)) 








z = 0.1 


z = 0.5 


z = 1 


z = 1.5 


z = 2 


El 


5.3 (0.2%) 


2.2 (13.6%) 


1.5 (34.9%) 


1.1 (45.0%) 


6.6 (2.5%) 


E2 


2.5 (2.4%) 


3.3 (7.4%) 


4.6 (4.9%) 


2.8 (16.4%) 


0.5 (80.8 %) 


B 


5.4 (0.3%) 


5.2 (3.2%) 


5.7 (4.0%) 


3.1 (12.2%) 


1.1 (50.8 %) 


H 


2.3 (2.5%) 


3.7 (6.0%) 


5.5 (4.3%) 


1.9 (26.9%) 


0.4 (90.8 %) 
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Figure 1. The AGN bolometric LF at 0.5 SC z 2: models vs. observations. The dashed black lines show the bolometric LF of H06, while the yellow shaded 
areas take account of the estimated errors of the fit. The dashed vertical lines show the minimum bolometric luminosity accessible to observations. Each set of 
plots, composed by four panels corresponding to different redshifts, refers to a different model, as indicated by the labels. 
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estimated the following % -like quantity: 



rmodel / 



1 



* [log(w mod ei(AL/,z)) -log(n obs (ALi,z))] 



bin ; =1 



(4) 

where n mo( j e ](AL,-,z) and n b s (AL,-,z) represent the model and ob- 
served mean comoving number density of AGNs in the luminosity 
interval AL; at redshift z, a m odel and a a ^ s are the la logarithmic 
errors and the sum runs over the N^ m luminosity bins in the interval 
^obs ~ ^model' ^ e nave verified that all our results are not sensitive 
to the choice of the bin size. 

The values of E model (z), evaluated at all redshifts, are shown 
in Table Q] for all models explored. To compare these values with 
theoretical expectations, we use Monte Carlo techniques to com- 
pute distribution of £, f(E,z), expected when « mot jei(AL,,z) is a 
Gaussian random variable, normally distributed around n a ^ s (AL/.z) 
with variance lO^ ™* 1 " 1 " " 1 "'. This function is used to evaluate the 
cumulative probability of S by integrating /(E,z): 



model 



z = 



dE/(E,z) 



(5) 



which is defined in analogy to the x 2 -probability and represents 
the probability that a function that genuinely describes the dataset 
would give a value larger or equal to s model . The values of 
P(E model ) z) are listed (in parentheses) in Table Q] and indicated in 
each plot. 

The results confirm those of the M06 analysis, in the sense 
that all models, apart from El that significantly overpredicts the 
abundance of AGNs at z = 2, match the LFs in the range 1 < z ^ 2 
fairly well. The advantage of considering bolometric rather than 
B-band or hard X-ray LFs is apparent at lower redshifts where dis- 
crepancies between models and observations at the bright and faint 
ends of the luminosity functions are more significant here than in 
the M06 analysis. Indeed, all models overpredict the abundance of 
bright objects and underpredict the abundance of the faint ones at 
z = 0.5 and z = 1. 

In the LF bright end, the mismatch can be reduced by ad- 
vocating some physical mechanism, like inefficient cooling, that 
hampers mass accretion in large haloes. Our simple model E2, in 
which mass accretion is inhibited in haloes with masses larger than 
10 13 ' 5 M©, provides a better match to data, especially at z = 0.5, 
although the effect is less apparent here than in M06 which consid- 
ered the optical B-band LF. The overabundance of bright AGNs is 
also alleviated in model H since the variable Eddington ratio guar- 
antees that a BH hosted in the largest halos accretes most of the 
time at a sub-Eddington rate, resulting in a fainter AGN. 

In all models, but El, the LF faint end is biased low. The effect 
is systematic and, in the luminosity range accessible to observa- 
tions, it does not depend on luminosity. Discrepancies grow larger 
when extrapolating the comparison below to objects fainter than 
^bof ' below which the LF predicted by most semi-analytic mod- 
els turns-over while the model LF of H06 is fitted by a power-law. 
The power-law behaviour is, however, recovered by model H, that 
assumes a time-dependent Eddington ratio. 



3.2 The hard X-ray luminosity function at z ~ 0. 1 

To understand whether the under-abundance of faint AGNs pre- 
dicted by most semi-analytic models is real or a mere artifact re- 
sulting from having extrapolated the power-law behaviour of the 
bolometric LF of H06 below L^Jjj 1 requires probing the AGN LF to 
lower luminosities, which is only possible in the nearby universe. 
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Figure 2. The AGN bolometric LF at z = 0. 1 : models vs. observations. S06 
and B06 LFs are represented by filled and open dots, respectively. Vertical 
error bars represent 10 uncertainties while horizontal bars indicate the size 
of the luminosity bins. Each plot refers to a different model, as indicated by 
the labels. 



In this section we do not compare the model LFs with the 
bolometric LF at z ~ 0. Instead, we apply the inverse bolometric 
conversion of H06 to compare model predictions with the LFs of 
S06 and B06 at z ~ 0. 1 in the [2 - lOkeV] band. The rationale be- 
hind this choice is as follow. First of all, these two datasets, espe- 
cially the B06 one, include objects that were not considered in the 
H06 analysis. Secondly, selection in the hard X-ray allows to in- 
clude obscured AGNs which make bolometric corrections less se- 
vere in this band. Third, the two samples have rather similar com- 
position as 90% of the objects are Seyfert galaxies. As a result, 
comparing model with S06 and B06 LFs allows to maximize the 
number of nearby, homogeneous objects, while reducing uncertain- 
ties in model bolometric corrections. 

Model vs. data comparisons are shown in Fig. [2] where the 
S06 and B06 LFs are represented by filled and open dots, respec- 
tively. Vertical errorbars represent la uncertainties, while horizon- 
tal bars indicate the size of the luminosity bins. The AGN lu- 
minosity in the B06 sample are measured in the [20 — 40keV] 
band and transformed in th e [2 — lOkeV] band according to 
^[2-l0keV]/^[20-40keV] = 2.3 teeckmann et alj200r3) . 

The shaded areas show the model LFs at z — 0. 1 together with 
their lo uncertainties expressed in the [2 — lOkeV] band by using 
the bolometric correction of H06 
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with sci = 10.83, ki = 0.28, c 2 = 6.08 and k 2 = -0.02. The bolo- 
metric luminosities are indicated on the X-axis in the upper part of 
the plot. In order to correct for the extinction in the X-ray bands, 
we have used the following equation, also provided by H06: 
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where = 1.243, P = 0.066 andL(,ol[^[2-l0keV]] ' s me bolometric 
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Figure 3. The AGN bias function at z < 2: models vs. obse rvations. The 
solid black points show the bias o f P orciani &Norberd i2006). The shaded 
areas show the bias predicted by our four models. The four dashed lines, 
with their shaded areas, show our model predictions with their la uncer- 
tainties. 



ISheth et aflfeOOlh , b(M^ a i ,z), with the mass function of the haloes 
hosting AGNs with luminosities larger than the thresholds of the 
observations, 1 !' (M, z) '■ 

f MMialo, z)V{M hllo {L B > L minjB ),z)dM hal0 
Mi) = - » ( 8 ) 

/ ^(Mhalo (L B > Anin^) , z)dM h . dio 
Jo 

where the minimum B luminosity at the 5 redshifts explored 
is L minM /L Q . B (z) = {1.42 • 10 n ,3.58 ■ 10 n ,3.92 • 10 n ,6.82 ■ 
10 u ,6.82- 10 11 , 1.08- 10 12 }. The difference with respect to the 
M06 analysis is that here we have used the bolometric correction of 
H06 to convert our bolometric luminosities to B-band ones, i.e. us- 
ing the Eq. © with c\ = 6.25, k\ = -0.37, c 2 = 9 and k 2 = -0.012, 
and the Eq. ifTJ with = 0.26, (3 = 0.082. We have compared 
model predictions with the most rec ent observational determina tion 
of the biasing function at ^ z ^ 2 flPorciani & Norberdl2006i) . es- 
timated in the B-band. 

The result is shown in Fig. [3] where the points represent the 
observed B-band AGN biasing function and the shaded areas show 
the lo" uncertainty interval around model predictions. As in the 
M06 analysis, the large model uncertainties do not allow us to place 
strong constraint based on the AGN clustering. Indeed, all our mod- 
els are in acceptable agreement with the data, suggesting, however, 
that possible disagreements may become significant at z > 2. 



luminosity correspondent to £<[2-lokeV]> as gi ven by the bolometric 
corrections of Eq. 

The comparison between model and data confirm our previous 
extrapolation, since the observed number density of faint AGNs 
with Lx = 10 42 — 10 43 erg/s is significantly larger than that pre- 
dicted by all models, as indicated by the sudden drop in the values 
of the p(s model ) at z = 0. 1. This is due to the fact that, for a given 
value of 2 model 5 the /(E) distribution at redshifts ^ 0.5 is more 
positively skewed than at z = 0.1, as we have verified. Discrepan- 
cies are larger for models El and B, while models E2 and H provide 
a better match to the faint end of the local LF. The sharp downturn 
in the El and B models is a robust feature since the characteristic 
mass of halos populating the faint luminosity bins (~ 10 Mq) 
is well above the mass resolution limit in our merger trees. 

Note that the largest discrepancies are found in the faintest lu- 
minosity bin which can only be pro bed by the B06 sample. With 
this respect, it is worth noticing that lSazonov & Revnivtsevl (2004) 
have used yet another dataset of hard X-ray selected AGNs to es- 
timate the AGN LF down to Lp_20keV] ~ 10 4I erg/s. Their LF is 
consistent with t hose of S06 and B06 down to the faintest objects. 
The sample of lSazonov & Revnivtsevl d2004h consists of 95 AGNs 
in the [3 — 20keV] interval at high galactic latitude serendipitously 
detected in the RXTE slew survey. However, only 60% of these 
AGNs are classified as Seyfert galaxies, many of which also belong 
to the S06 sample. Since in this work we prefer to deal with a ho- 
moge neous sample of local AGNs , we have decided not to consider 
the Sazonov & Revnivtsevl J2004l) LF in our quantitative analysis. 

3.3 The biasing function 

For the sake of completeness, we follow the M06 analysis and 
quantify the clustering of our model AGNs at z ^ 2 through their 
biasing function, b e ff(z). The latter has been computed by weight- 
ing the analytic biasing function of the DM haloes provided by 



4 DISCUSSION AND CONCLUSIONS 

In this paper we tested the validity of the assumption that the evolu- 
tion of AGNs is simply related to the cosmological merging history 
of DM haloes. To do that, we have compared the predictions of hi- 
erarchical semi-analytic models with the most recent determination 
of the AGN LF in the hard X-ray band and their biasing function 
in the B-band at z ^ 2. Our main results can be summarized as 
follows. 

(i) We confirm the success of simple semi-analytic models in 
reproducing both the AGN bolometric LF at 1 < z ^ 2, i.e. around 
the peak of activity, and their clustering, quantified by the biasing 
function, at z ^ 2. 

(ii) As pointed out by several previous analyses, problems oc- 
cur at lower redshifts, where hierarchical models systematically 
overestimate the number density of bright AGNs and underestimate 
the faint ones. 

(iii) Comparing bolometric LFs rather than the optical or hard 
X-rays ones allows to spot significant discrepancies already at mod- 
erate redshifts z ~ 1, i.e. earlier than what was found in previous 
analyses (e.g. M06). 

(iv) The predicted number density of bright AGNs can be re- 
duced not only by advocating inefficient cooling within massive 
haloes, as in model E2, but also by accounting for feedback mech- 
anisms, as we did in model H. 

(v) The underestimate of faint AGNs looks like a more se- 
rious problem that we have tried to tackle by assuming a time- 
dependent Eddington ratio, as sugge sted by the outc ome of the 
hy drodynamical s i mulati ons by iHopkins et alj {2005). As shown 
bv lVolonteri et alj J2006t) . implementing this prescription within a 
semi-analytic framework, as we did in model H, proved to be suc- 
cessful in reproducing the r edshift distribution of the faint X-ray 
counts JVolonteri e t al. 2006). In this work, we extended the analy- 
sis o f IVolonteri etail J2006h by comparing model predictions with 
the most recent determinations of the local AGN LF by S06 and 



Hierarchical evolution ofAGNs 7 



B06 in the hard X-ray band, to include absorbed AGNs and mini- 
mize the impact of bolometric corrections. This is a very demand- 
ing test for semi-analytic models, which constitutes the main focus 
of this paper since, as we have pointed out, the mismatch in the 
number density of faint AGNs grows larger when decreasing the 
redshift. We found that the two most successful models E2 and H 
are in acceptable agreement with the data at z > 0.5, but struggle 
to match the correct number density of faint X-ray sources in the 
nearby universe. 

Mo del H, based on the results of hydrodynamical simula- 
tions of Hop kins et all J2005h within a pure merger driven sce- 
nario, seems unable to account for local faint AGNs. If the ac- 
cretion efficiency were much lower, the lifetime of faint AGNs 
would increase proportionally and help alle viate the discrepancy . 
However, the Eddington factors derived from Hopkins et al. (2005) 
light curve are well below /gdd = 0.1 only when a galaxy hosts a 
black hole with an initial mass anomalously smaller than that pre- 
dicted by the Mbh — a g correlation. This is evident in Figure [4] 
our models assume that accretion processes bring the black holes 
onto the Mbh — <5g relation and the accretion efficiency is for most 
systems above /gdd = 0.1. This can be understood using a very 
simple argument. Let us assume that (i) quiescent black holes sit 
on the Mbh — <5g relation, as observed in the nearby galaxy where 
the Mbh — a g relation was indeed derived. This is therefore a safe 
assumption in the local Universe, (ii) Accretion is triggered only 
by major mergers, that is mergers betw een similar size galaxies, 
with a mass-ratio larger than at least 0.1 (Cox 2004). And, (iii) an 
accretion episode grows black holes until they reach the Mbh — Gg 
relation for the newly formed galaxy, due to feedback effect. Stud- 
ies of local samples of AGN seem indeed to confirm that typically 
AGN masses scale with the Mbh — °g relation (e.g., iFerrarese et al] 
1200 ll ; iGreene & Hoi [2006). Within these simple but sensible as- 
sumptions, the accretion efficiency is bound to be high, as can be 
easily estimated by equation [3] If we consider, for example, a ma- 
jor merger of a Milky-Way sized galaxy, the Eddington factor of 
the black holes remains /^dd < 0.1 for only about 10 6 yr. 

The inadequacy of the pure merger driven scenario becomes 
more evident when considering the observational constraints on the 
Eddington factor of Seyfert galaxies, which con stitute about 90% 
of the local AGN population. IWoo & U rrv ( 2002) analyze a sample 
of 234 AGNs at 0.001 < z < 1, composed, at z < 0.1 mainly by 
Seyfert galaxies. They find a large scatter (2 orders of magnitude) 
in the Eddington fac tor at both fixed luminosity and fixed BH mass. 
IWoo & Urry] j2002l) do not find any trends of the Eddington factor 
with either luminosity, mass or redshift, which cannot be explained 
by selection effects. 

It turns out that the S06 and B06 catalogues are largely com- 
posed by Seyfert galaxies that constitute respectively 94% and 88% 
of the total galaxy host populations. Only a small fraction of these 
local Seyfert galaxies have disturbed morphology, and thus did 
not experience any recent merging event. Indeed, only 4% of the 
sources in the S06 catalog are hosted in galaxies that show ev- 
idences of recent interactions. The AGNs in the B06 catalogue 
are typically found at low galactic latitudes which hamper a sys- 
tematic analysis of their host galaxy morphology. Yet, the similar 
galaxy composition of the two catalogues suggests that also B06 
AGNs preferentially populate quiescent environments. Based on 
this observational evidence, it may be suggested that galaxy merg- 
ers might not constitute the only trigger to AGN activity. 

To decide whether this is indeed a viable hypothesis, it is 
worth reviewing the observational evidences of local AGN sam- 
ples. Bright, low-redshift quasars and ultra luminous infrared 




8 10 12 14 



log(L bo ,/L lolE ) 

Figure 4. The mean Eddington ratio in function of the AGN luminosity, at 
z = 0.1, for the model H. The coloured area represents the Io~ uncertenties. 

galaxies, that are generally regarded as hosting obscured AGNs, 
are often found in merging systems (see, e.g. ISanders & M irabel 
ll996tlCanalizo & StocktorfeOOlHCapetti & BalmaverdJ 20061. and 
references therein) which indicates a possible connection between 
mergers of gas-rich galaxy and AGN activity. On the contrary, 
as we have seen, fainter AGNs typically reside in quiescent, 
non-i nteracting galaxies (e.g. lKauffmann etal ] |2003l ; lGroginetal.l 
2005, and references therein). However, this alone does not guaran- 
tee that an alternative AGN triggering mechanism is at work, as this 
observational evidence can still fit into a merger-driven scenario. In 
fact, the brightest among these objects could be the relics of a pre- 
vious bright quasar epoch in a spheroid-forming merger (see, e.g. 
Hopkins et al. 2006, and references therein), while the fainter ones 
would consist of AGNs hosted in "dead" elliptical galaxies fueled 
via accr etion of hot spheroid gas and steady mass loss from stars 
(see, e.g. lciotti & Ostrikej200ll ; ISazonov et al]2005l : lcroton et al] 
2006), an accretion mode which cannot dominate the BH growth. 

The merger-driven scenario, however, proved to be inade- 
quate in accounting for the relatively high accretion rate AGNs 
observed at low redshifts in undisturbed, late-t ype, star-forming 
galax ies with low mass (< 10 7 M Q ) BHs (e.g. Kauffmann et al. 
120031) . Indeed alternative mechanisms, not included in our sim- 
ple models, have been suggested to trigger the mass accretion in 
these objects. For instance, it has been proposed that a signif- 
icant contribution to the faint AGN mass accretion could come 
from the materia l liberated by the tidal d isruption of stars by 
the central BHs (Milosavlievic et al. 2006), or by the mass of 
the stars captured by the BH dis ks and eventually dissolved 
dMiralda-Escude & Kollmeiej l2005h . Other studies have consid- 
ered the stochas tic accretion of molecular clouds in quiescen t 
systems (see e.g. Hopk ins & Hemauisdl2"006r . ICroton et al] |2006). 
Moreover, it was suggested that also disk instability could trig- 
ger mass accretion, contributing t o increase even more the num- 
ber de nsity of faint AGNs (see e.g. lCroton et al] |2006; Bow er et"al] 
2006). Finally, a better treatment of mergers between haloes with 
low mass ratio may also contribute to s olve these problems (see e.g. 
iMalbon et alj|2006r . lCroton et al]200fj) . 
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Whether including these alternative trigger mechanisms in our 
simple merger-driven scenario can help in reconciling model pre- 
dictions with observations at z ~ is a question that deserves fur- 
ther investigation. 
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